Dynamical chiral symmetry breaking in SU(N C ) gauge theories with large number of 

fermion flavors 
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In this paper we examine a phase transition in SU (N c ) gauge theories governed by the existence of 
an infrared fixed point of the renormalization group /3 function. The nonlinear integral Schwinger- 
Dyson equation for a mass function of massless fermions is solved numerically using the exact 
expression of the running coupling in two-loop approximation for an SU(3) gauge theory. Based 
on the obtained solution of the Schwinger-Dyson equation, the value of the chiral condensate, (qq) , 
and the decay constant, /„■, of bound states (mesons) are calculated for several values of fermion 
flavors Nf . We show that this kind of phase transition is a transition of finite order. 
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I. INTRODUCTION 

Gauge field theories with a large number of massless 
fermions are becoming an attractive topic for theoretical 
research. A recent discovery of a phase transition in such 
theories [l|, @ has led to additional interest. Depend- 
ing on the number of fermions, there are two possible 
phases. The first phase is a phase with broken chiral sym- 
metry and confinement which occurs when the number 
of fermions is less than some critical value (Nf < Nj r ), 
where Nf r is a critical value of the number of fermions. 
The second phase, occurs when (Nf > N^ r ), is a phase 
with strict chiral symmetry and absence of confinement. 
This type of phase transition occurs for example in min- 
imal supersymmetrical QCD [3| . The dynamics of these 
two phases is well studied in the approximation discussed 
bellow. 

The reason for the phase transition is the existence 
of an infrared fixed point (IFP), with coupling constant 
a* = a(0), in two loop approximation for the renormal- 
ization group (3 function. Also the mathematical aspect 
of the phase transition becomes clear from the analysis 
of the integral Schwinger-Dyson equation (SDE) for the 
mass function. Assigning to the running coupling the 
constant value at the IFP, a(Q 2 ) = a*, the SDE turns 
out to be an equation for eigenvalues and has only trivial 
solution with sufficiently small a*. In contrast, if a* is 
larger than a critical value a c = tt/(3C2(F)), then the 
SDE has nontrivial solutions. The critical value N^ r is 
then determined from the relation a* = a*(Nf,N c ), at 
constant N c . However the situation changes when we use 
the running coupling. The kernel of the integral SDE, 
like the running coupling, is a function of the number 
of fermions and colors. And it is unclear what should 



be considered as a critical value of a*. The most obvi- 
ous solution would be the introduction of a certain in- 
tegral characteristic for the SDE (as has been done for 
the symmetric kernel in the theory of integral equations). 
We could then consider number of fermions as the free 
parameter, which would determine the value of a*. 

If a chiral symmetry is broken then there exist boson 
degrees of freedom, which of course arise as Goldstone 
bosons. We examine the chiral phase transition by study- 
ing order parameters like the quark condensate, (qq), and 
the decay constant, f„, of bound states. Examining these 
quantities near the point of phase transition may shed 
light on the nature of the transition. In this paper we 
study numerical solution of the SDE for a mass function 
using an exact expression for the running coupling. The 
existence of both trivial and nontrivial solutions of the 
SDE confirms that the phase transition takes place in an 
SU(3) gauge theory. 

In section [TT] we discuss the equation for the running 
coupling in two loop approximation. In section IIIII we 
briefly go over the conservation of Ward-Takahashi iden- 
tities and set up the SDE in the local gauge. Section ITVl 
is devoted to the numerical calculations and discussion. 



II. SOME PROPERTIES OF THE GAUGE 
THEORIES WITH THE INFRARED FIXED 
POINT 

Let's start from the Lagrangian of the SU(N C ) gauge 
field theories. It appears as follows: 
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where tpk is a four component spinor of flavor k, = 
<9 M — igA^T a , T a are the generators of the gauge group 
and g is the coupling constant. This Lagrangian is clearly 
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invariant under the global symmetry group SU(Nf)i, x 
SU(Nf)ji x U{1)l,r because all fermions are massless. 
Nonetheless, this symmetry may be broken to diagonal 
subgroup SU(N f ) L+R x U(1) L , R . 

The next step is an analysis of the equation for the 
running coupling in two loop approximation (the first 
two coefficients are independent of the renormalization 
scheme, the higher-order coefficients are scheme depen- 
dent). It takes the form: 
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where a — g 2 /4tt and according to Ref. Q coefficients b 
and c look as follows: 



b=—(llC 2 (A)-4T f N f ), 



(11.3) 
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The theory is asymptotically free if b > 0. The IFP exists 
if c < 0, which for SU(3) takes place when Nf > 8. The 
running coupling at the IFP takes the value a* = —b/c. 
The fixed point coupling a* can be made sufficiently 
small to perform a calculation by perturbation theory. 
Certainly we assume that < g < 1 in the initial La- 
grangian (|II.1[) . which restricts the variation interval for 
Nf . For an SU (3) gauge theory this variation interval is 
10 Nf < 33/2 assuming an asymptotic freedom of the 
theory. 

The equation (|II.2[) can be integrated as follows: 
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where we have introduced the scale 



A = v 



a 8 (v 2 ) 



c/b 

c ) exp 



Q 2 
6 In -pr, 

A 2 



ba(v 2 



(II.5) 



(II.6) 



which has the same physical sense as the dimensional 
Aqqd parameter in ordinary QCD. For further calcula- 
tions we consider A to be independent of Nf, for a fixed 
value of N c . This is a good approximation when we are 
taking into consideration the small variation interval of 
Nf. 

Equation (|II.5jl is a transcendental equation for a(Q 2 ). 
It can be solved analytically using complex Lambert 
Wk(z) function Lambert's function satisfies the 

transcendental equation W^(z) exp(Wfc(z)) = z, where 
k = 0, ±1, ±2, .... We note that there are only two pos- 
sible real solutions k = 0,-1 which take place if the ar- 
gument z — 1/e. The Lambert W function has simple 
asymptotics: Wo(z — > 0) ~ z — z 2 , W-\ diverges when 
z — > 0, and Wk{z — > exo) ~ lnz + 2i:ik — ln(lnz + 2itik). 



Then we have: 
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The case i = — 1 corresponds to the running coupling in 
ordinary QCD in two loop approximation (in the pertur- 
bative regime). Further, only the case c < will be used 
for calculations. Based on the properties of Lambert's 
function, it is possible to obtain the asymptotic form of 
(HIT]): 



a(Q -> 0) = a* 1 + 
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It is important to note that as the number of fermions 
decreases, the value of asymptote (|II.8j) increases, and 
the value of asymptote (|II.9[) decreases. 



III. SCHWINGER-DYSON EQUATION FOR 
THE MASS FUNCTION 



Since the initial Lagrangian has only massless 
fermions, we must require conservation of the vector and 
axial vector Ward-Takahashi (WT) identities. It is essen- 
tial for the conservation of the axial WT identity that the 
running coupling must depend on the same momentum 
as the gluon propag ator 181. The use of simple Landau 
gauge and other local covariant gauges violates the vec- 
tor WT identity when the ladder approximation with a 
bare vertex is studied, and they are not suitable for this 
purpose. This problem can be solved by using a nonlocal 
gauge which depends on the momentum: 
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This nonlocal gauge was proposed by T.Kugo and 
M.Mitchard f9], where r)(p): 
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It is clear that function (|III.2j) coincides with the Lan- 
dau gauge at small and large momenta, i.e., 77(0) = 
1, 77(00) = 1. 

The nonlocal Kugo gauge allows us to write the SDE 
for the dressing fermion propagator S(p) — i/(A(p 2 )p — 
B(p 2 )). In the ladder approximation it is given by: 
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where we use a bare vertex ig^^T a . The another advan- 
tage is that in the Kugo gauge the fermion wave function 
renormalisation constant A(p 2 ) is equal to one. Using 



(IIII.3I) we may retrieve the SDE for the dynamical mass 
function in Euclidean momentum space: 
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However, the obtained equation is rather complicated 
and can not be solved analytically without certain as- 
sumptions and approximations. One of these approxi- 
mations is to use the constant value at the IFP for the 
running coupling, i.e., a(Q 2 ) — a*. It may be used in 
the region of momentum where the running coupling is 
slowly changing. Or in other words for a* — > a c from 
above and Nf — > N c from below. Expression (|III.2p in 
this approximation is equivalent to the Landau gauge and 
therefore (|III.4[) becomes: 
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Equation (jHI.5j) can be integrated over angular part: 
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In general this type of angular integration can not be 
performed exactly in the more complicated cases. Sub- 
stituting (|III.6|) into piI.Sp . we obtain the final integral 
equation for the mass function: 
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were we have introduced the cutoff parameter A* . By us- 
ing simple differentiation with respect to p 2 , this integral 
equation converts to a differential equation: 



p 2 B'\p 2 ) + 2B'(p 2 ) 



3a*C 2 (F) B(p 2 ) 



An 



B 2 {p 2 ) 



with two boundary conditions 
dp 2 B(p 2 ) 



= 0, 
(III. 
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The solution of this equation can be found in terms of a 
hypergeometric function which has the following asymp- 
totic form at small momentum: 



B(0) rj A„exp 



C 



\/a*/a c 



The origin of the critical constant a c now becomes fairly 
clear. For recent reviews of the SDE and their application 
see for example Refs. @, Q. 



IV. NUMERICAL CALCULATIONS AND 
RESULTS 



In this section we discuss the method which was used 
to solve the SDE (|III.4[) . First of all we note that it is im- 
possible to perform the angular integration analytically 
in this case. To solve the nonlinear SDE, we use a sim- 
ple quadrature method. All calculations were performed 
using Mathematica software and consisted of the follow- 
ing steps. We set up a square lattice (pi, kj) where both 
(III 10) Pi an d % P ass the number of discrete values from lower 
boundary qo to upper boundary A* and then carry out 
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FIG. 1: Numerical solution of the Schwinger-Dyson equation. 
Here B/A is dimensionless dynamical mass function 



TABLE I: Numerical values of the physical quantities ob- 
tained from the numerical solution of the SDE 



N f 


B(0)/A 


-(qq)/A 3 ■ 10- 4 


/tt/A 


10 


0.0951 


5.3224 


0.02791 


10.3 


0.0540 


2.2542 


0.01845 


10.5 


0.0336 


1.0331 


0.0126 


10.7 


0.0186 


0.3686 


0.0077 


11 


0.0055 


0.0473 


0.0028 



numerical integration using these lattice sites. Next, we 
replace the integral in (|III.4[) by quadrature sum of the 
rectangles. As a result we obtain a system of nonlinear 
equations, where unknown variables serve as the values 
of the unknown mass function at the lattice sites: 



B(pi) 



kj —q 



A, (IV.l) 



where A is the lattice step. 

Since the unknown function is smooth, it is possible 
to use a quadrature formula of rectangles. Also, the 
same results are reproduced when we use a more pre- 
cise quadrature trapezoid formula. The newly obtained 
system of nonlinear equations was solved by an iteration 
method. We have calculated the mass function for a se- 
ries of values of Nf near Nf . 

In Fig[T]we illustrate the dependence B/A on p/A for 
various numbers of fermions. The behavior of the mass 
function is as expected: it differs from zero at small 
momentum and then go smoothly to zero at large mo- 
mentum. Also the quantity B /A is obviously dimension- 
less and cannot depend on the transmutation parame- 
ter A. The mass function B increases when the number 
of fermions decreases. There is a similar tendency in 
(IIII.10|) . The cuttoff A* is approximately equal to 0.3A. 
In this area the approximation a(Q 2 ) ~ a* is useful and 



it is interesting to compare B(Q) with (|III.10|) . We will 
do this below. 

For better understanding of the physical nature of the 
phase transition, it is also useful to calculate other phys- 
ical quantities. One of them is the value of the vacuum 
condensate, which can be easily found: 
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and the value of meson decay constant (it is described by 
well known Pagels-Stokar formula) 
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B(p) 



P 2 B(p) 



dp 2 
(IV.3) 

These values have also been calculated using quadrature 
formulas and are illustrated in Table [U 

Let us analyze the quantities -B(O), (qq), f v in more de- 
tail. These quantities are continuous functions near the 
critical point. This fact confirms that the phase tran- 
sition is a transition of second or higher order, possibly 
of infinite order. For this reason, the quark condensate, 
decay constant and B(0) may be fitted by polynomial 
functions near the critical point. 
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If we find that the critical exponents a, (3, p are small 
real numbers, it will confirm that the phase transition is 
of finite order. Note, that Eg. pil,10[) describes a phase 
transition of infinite order as oppose to finite. The least 
squares fitting of the curves gives: 



B(0): 


a 


= 2.45, 


Nf 


= 11.43; 
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= 2.93, 


Nf 


= 11.21; 
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= 1.69, 


Nf 


= 11.33. 


(IV.9) 



We indeed find that the critical exponents are small 
positive numbers. The value of the critical number of 
fermions is described by the well known formula: 



Nf = N c 



(WON 2 
\25Nf- 



66 
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obtained from the condition a* and a c [l| . In case N c 



3, the critical number of fermions is N c t 



11.9. 



V. CONCLUDING REMARKS 

In this paper we have shown numerically that the SDE 
for a fermion propagator with an exact expression for the 
running coupling has nontrivial solution B{p). Based on 
the obtained solution, the value of chiral condensate and 
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the decay constant of pseudoscalar bosons were calcu- 
lated. These physical quantities are continuous functions 
near the critical point. Detailed analysis of B(0) near 
the critical point shows that the phase transition is a 
transition of finite order. 
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